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Abstract: This paper considers Darboux transformations of a bispec- 
tral operator which preserve its bispectrality. A sufficient condition 
for this to occur is given, and applied to the case of generalized 
Airy operators of arbitrary order r > 1. As a result, the bispec- 
trality of a large family of algebras of rank r is demonstrated. An 
involution on these algebras is exhibited which exchanges the role 
of spatial and spectral parameters, generalizing Wilson's rank one 
bispectral involution. Spectral geometry and the relationship to the 
Sato grassmannian are discussed. 

Introduction 

An ordinary differential operator, L(x, d x ), is said to be bispectral j$J if it 
has an eigenfunction ip(x, z) satisfying a pair of eigenvalue equations 

L(x, d)if>(x, z) = f(z)if>(x, z) Q(z, d z )ip(x, z) = g{x)i>(x, z) (1) 

for non-constant functions / and g. This property was investigated by 
Duistermaat and Grunbaum ||/, who identified all bispectral Schrodinger 
operators: L = d 2 + V(x). In j7^/, Wilson considered the more general 
question of classifying bispectral commutative rings of ordinary differential 
operators. That is, he considered rings £ and Q of operators in the variable 
x and z respectively such that there exists a function if}(x,z) for which 
Equation (|I|) holds for all L e £ and Q G Q. 
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Following j^j, one defines the rank of a commutative ring of ordinary 
differential operators to be the greatest common divisor of the orders of its 
elements. Wilson classified all maximal bispectral rings of rank one: 

Theorem 1.1: [Wilson, Let 1Z be a maximal rank one commutative 
algebra of ordinary differential operators. Then 1Z is bispectral if and only 
if its spectral curve is a rational curve with no singularities other than 
cusps. 

This result was proved, in part, through the introduction of an invo- 
lution on a subset Gr ad of the grassmannian Gr l generally used in the 
investigation of the KP hierarchy /|73[/. By the work of Sato, Segal, Wil- 
son and others \ 7^] one knows how to associate to a point W G Gr 1 a 



commutative ring of operators and the Baker function tpw(x, z ) which is a 
common eigenfunction of the operators in the ring. Then Gr ad is the sub- 
space of Gr 1 consisting of suitably normalized points whose corresponding 
ring is the affine coordinate ring of a rational curve with cusps, minus a 
point at oo. Wilson showed that Gr ad is equipped with an involution, 
j3 : Gr ad — > Gr ad , with the property that 

if>p(W)(x, z) = ip w {z,x) (2) 

and consequently that the corresponding rings are bispectral. (Essentially, 
this is the involution which exchanges the rings £ and Q.) 

The classification of bispectral algebras in higher rank appears to be 
more subtle |^/. There is, however, a well-known method for construct- 
ing new commutative algebras of differential operators from a given one: 
Darboux transformation. (See, for example, Jjldj , \FJj , \H\] for previous results 
relating these transformations to the bispectral problem.) Under suitable 
hypotheses, described below, Darboux transformations preserve bispectral- 
ity In this paper we will explore this method in the case of generalized 
Airy operators. 

The general approach is the following. Let T> denote the ring of ordinary 
differential operators with coefficients in C(x). Given a function ip(x,z), 
consider the equation 

T(x,d)^ = T\z,d z )^ , (3) 
for T and T b in T>. Define 

A^CV (4) 

to be the set of operators T for which there exists such that Equation (^) 
holds. Note that is a subalgebra of V. If i/j is sufficiently general, in 
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particular if -0 is a bispectral eigenfunction satisfying equation ([!]), then 
the map 

T ^ T b (5) 
is an anti-isomorphism from A^ to A^, where 

i]){x, z) = if)(z, x) . (6) 

Given a pair of operators T and L, define an algebra of polynomials 

Rt,l := { p{z) e C[z] | Tp(L) COT}. (7) 

Now assume Equation ([!]). For reasons explained in /[73|/, it suffices to 
assume that L and Q are normalized, by which we mean that the top two 
coefficients are, respectively, constant and zero. It then follows that L and 
Q have coefficients in the field of rational functions C(x) (resp. z), and 
that f(z) and g(x) are polynomials. For T G A<p, let a(x) (resp. a b (z)) be 
the leading coefficient of T (resp. T b (z,d z )) and 

f:=^-T f b :=^-T\z,d z ). (8) 

Then we immediately have the following proposition. 
Proposition 1.2: Given pi G Rt,l and P2 G Rt^q, ^ 

P 1 :=fp 1 (L)f- 1 (9) 
P2 := f^iQ)^ 1 (10) 

and let 

<K?,*)-= f \ ( M (x,z)) . (11) 
a{x)a {z) 

Then Pi and P2 are normalized ordinary differential operators, satisfying 

Pi<f> = pi{f (*))<!> (12) 
P2(z,d z )<P = p 2 (g(x))(P . (13) 

Corollary 1.3: Given a bispectral triple (L,Q,ip) satisfying Equation^, then 
for T G Aip such that the rings Rt,l and Rt\q are non-trivial, the ring 

TZ T , L := {fpi^f- 1 I p G R t ,l } (14) 
is an algebra of normalized bispectral operators, isomorphic to Rt,l- 
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Thus one has a strategy for obtaining new bispectral algebras from a 
given bispectral operator, namely to find T G A^ such that the rings Rt,l 
and Rp> q are non-trivial. Each such T then gives a bispectral Darboux 
transformation of the bispectral operator p(L) for p G Rt,l- 

For instance, Wilson's result, though not stated as such, is essentially 
the following. Taking ip(x, z) = e xz , A^ is the Weyl algebra of differential 
operators with polynomial coefficients. Indeed, in that case the map T i— > 
T b is the standard antiautomorphism of the Weyl algebra exchanging d 
and x. 

roposition 1.4: Let X 1 , X n be (not necessarily distinct) complex numbers 
and let pi(x) , . . . , p n (x) be polynomials. Let (pi = pi(x)e XiX . Let K be the 
operator 

K(u)x=e^*\Wr{fa,...,<l> n ,u)\ , (15) 

where \Wr((f)i, ...,(f> n ,-)\ denotes the Wronskian operator. (Clearly K has 
polynomial coefficients.) Then 

1. The rings Hjtd an d ^-R b d are both nontrivial. 

2. Every maximal, normalized rank one bispectral algebra is obtained in 
this way. 

The purpose of the present paper is to extend part of the proposition 
above to the higher rank case. What we are able to prove at the present 
time is that if we replace d by a generalized Airy operator 

L := d r - a r _ 2 d r - 2 a x d-x (16) 

for r > 1 and a$ G C, then the analogue of part one of the proposition 
continues to hold, at least when Ai,...,A n are distinct and pi,...,p n are 
polynomials of degree one. This analogue is the following. The operator 
L is clearly bispectral. Indeed, if we define 

f(x,z) :=f(x + z) , (17) 

then any / G ker(L ) satisfies the equations 

Lq{x, d)f(x, z) = zf(x, z) L (z, d z )f(x, z) = xf(x, z). (18) 

Note that for nonzero / belonging to ker(Lo), Af contains d, since d(f) = 

d z (f), and it contains x by equation (|T8|) . Thus At is again the Weyl 
algebra and so L$ determines an involution of the Weyl algebra. Explicitly, 
this involution is the unique anti-automorphism given by 

a> = L 9 b = d . (19) 
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Throughout the remainder of the paper we will consider L fixed and define 



the involution b by Equation 19 



The main result of the paper, to be developed and proved below is 

Theorem 1.5: Let fi(x), f r {x) be a basis for the kernel of Lq. Define K 
to be the Wronskian operator of the rn functions 

<ki - ■ (20) 

For Ai,...,A n distinct and pi, ...,p n polynomials of degree one, the rings 
^■k,l an d ^Ki> l are b°th nontrivial and therefore bispectral. 

One of our purposes here is to draw attention to the rather interesting 
subset of the Weyl algebra consisting of those elements T for which both 
1Zt,l and TZtKLq are nontrivial. We conjecture that it corresponds by the 
construction above to the set of homogeneous finite dimensional subspaces 
of the space of finitely supported distributions in the complex plane, as in 
Wilson's case. 

By a finitely supported distribution we mean a finite linear combination 
of operators of the form d~\odi, A G C, where 5\ is the 5-function evaluating 
its argument at z — A. We call such a distribution homogeneous if it is 
supported at one point, and we call a space of distributions homogeneous 
if it has a homogeneous basis. 

Let S be the space of finitely supported distributions in the z-plane, and 
let Grh(S) denote the set of finite dimensional homogeneous subspaces of S. 
Grh{S) appears to be a rather interesting space from an algebro-geometric 
point of view. Proposition |1.4| implies that Wilson's bispectral involution 



takes place on a certain quotient of Gr^S). Notice that if equation (0) 
holds, then L itself belongs to A^. Moreover, if T is replaced by Tq(L), 
where q(z) is any nonzero polynomial, then the rings Rt,l and R T \, L remain 
unchanged, as do their script counterparts. Now it is not hard to show that 
if C G Grh(S), and A is a complex number such that 5\ G" C then 

K C(B{ s x) =K c{z -x)o(d-\) , (21) 

where K now refers to the operator in Proposition TA. Thus if C G Grh{S) 



is a space of distributions containing a delta function, we can use the 
action of C[z] on Gr h (S) to replace point C with a point C' that contains 
no 5 functions and produces the same bispectral algebras. Let us call C 
minimal if it contains no 5 functions. Then every point C has a minimal 
representative in the above sense, so we can think of the minimal points 
as either a quotient or a subset of Grh(S). Proposition [O] implies that 
the rank one bispectral involution may be regarded as an involution on 
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the minimal points of Grh(S). Moreover, the involution takes place on 
finite dimensional pieces of Gr h (S), cut out by placing an upper bound on 
both the order K and the degree of its leading coefficient in Proposition 



L4j. We conjecture that precisely the same sort of phenomenon occurs in 
the generalized Airy case, though here we have only established this fact 
for elements of Grh(S) having a certain form. The natural approach to 
extending these higher rank involutions to all of Grh{S) seems to consist of 
first making a study of Grh{S) as an infinite dimensional algebraic variety 
and then applying a continuity argument. This we hope to do in a future 
work. 

1.6 Notation 

The notation ":=" will be used to indicate an initial definition of the ob- 
ject on the left hand side. Angular brackets (q) indicate the linear space 
spanned over C by the basis elements q. The determinant of a square 
matrix M is denoted \M\. The Wronskian matrix Wr(v) of an n vector 
v is the n x n matrix whose first row is v and so that each row is the 
derivative of the previous row. The symbols d t will be used to indicate 
and d is just d x . An ordinary differential operator L in the variable x is 
a polynomial in the symbol d with coefficients which are functions of the 
variable x. Often, a differential operator will be indicated by the notation 
L(x,d x ) to indicate an operator in the variable x and L(z,d z ) to indicate 
the same operator following the change of variables x — ► z. 

2 The Operators K and i? b 

First we prove a general lemma which will be used at several points through- 
out the paper. Note that our distributions are assumed to act in the z- 
variable, so that if -0 is a function of x and z, c(ifj) is a function of x. 

Lemma 2.1: Let L(x,d) be an rth order differential operator with coeffi- 
cients analytic in a neighborhood U C C. Let V C C be an open subset 
and let ipi(x, z), . . . , "0 r (x, z) be functions analytic inlA x V such that 

1. L(^j) = zipi for % = 1, . . .r. 

2. For all A G V, ipi(x, A), . . . , ip r (x, A) are linearly independent func- 
tions of x. 

Let OiU) be the space of analytic functions on U and S(V) C S be the 
space of distributions with support in V . Then the map 

S r {V) — > 0{U) 
( Cl ,...,c r ) i * E c *(^) ( 22 ) 
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is injective. 

Proof: Arguing by contradiction, assume Y,<k{i>i) — 0, with at least one 

i 

Ci non-zero. Let Ai, . . . , A n be the points at which the c's are supported. 
For j = 1, . . . ,n, let fj,j be the highest order derivative occurring among 
the c's at Xj. Then let m(z) be the polynomial 

n 

m(z):^(z-X 1 )^H(z-X i r i+1 . (23) 

i=2 

Let := q o m(z). Then = aid\ 1 with ai ^ if and only if q has order 
//i at Ai. In particular, c' l5 . . . , d r are not all zero. However, we have 

= m(L)(5>(V>i)) 

j 

= ^qom(L)(^) 
= ^ q o m(z)(ipi) 

i 

= J2 a ii>i( x , X i) ( 24 ) 

i 

This implies that the oj's are zero, which is a contradiction. □ 

Fix a basis < j < r} for ker Lo- Since the (r — l) s * coefficient of Lo 
is 0, the Wronskian determinant, \Wr(fi, . . . , f r )\ is a non-zero constant. 
Assume the basis has been chosen so that 

\Wr(f 1 ,...J r )\ = l. (25) 



Definition 2.2: Define D C Grh(S) to be set of all finite dimensional 
subspaces C C S having a basis of the form 

:= S Xt o (d z + 7j)|Aj distinct}. (26) 

Note: The elements of D should be thought of as identifying a singular 
rational spectral curve with simple cusps at the points Aj and a line bundle 
given by (71, ... , 7„) in the generalized Jacobian of that curve. 

Fix an element C E D. For convenience, we will order the basis for C, 
and define 



$ '■= (ci(/i),...,ci(/ r ),...,c n (/i),...,c n (/ r )) 

:= (M*),---,Mx)) (27) 
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where N := rn and f(x,z) = f(x + z). As an immediate corollary of 
Lemma |2.1| we have 

Proposition 2.3: The N functions, S; are linearly independent. 

The operator K in the introduction is defined only up to a constant, 
as it depends on a choice of basis. To be precise we will define 

Ku{x) = - ^ — \Wr{lu)\ , (28) 

where the sign will be specified in Theorem |2^4| below. 

It follows from Proposition |273| that |W / r(0)| ^ 0, so that K is indeed 
a differential operator of order N. We may thus consider the ordinary 
differential operator K 

Ku( X) = < 29 > 

\Wr(<f>i,...,<j> N )\ 

Note that this is the unique, monic operator of order N having the space 
(4> s ) as its kernel. 

The main result of this section is 

Theorem 2.4: 

1. K has polynomial coefficients. 

2. We can (and do) choose the sign of K so that the coefficient of d N 
is a monic polynomial of degree n. 

3. also has degree N . 

Before proving this theorem, it is convenient to modify the definition 
of the Wronskian matrix to take account of the present situation. Define 
a sequence of operators 

Pj = d l L k Q , where j = rk + l, < I < r. (30) 

Then Pj is monic of degree j. Given any vector of functions v = (vi(x), . . . , v, 
define the modified Wronskian matrix 



Wr{v!, ...,v m ):-- 



( * 
P x v 

V Pm-lV J 



(31) 



Since {Pj} is related to {<9 J } by a unipotent transformation, we clearly 
have 

\Wr(v u ...,v m )\ = \Wr(v u ...,v m )\. (32) 



Alex Kasman and Mitchell Rothstein 9 



Lemma 2.5: Given c = 5\ o (ad z + b) 6 S, 

Wr{c{h), c(/ r )) = (bl + aB(x + X))Q(x + A) 
where Q(x) := Wr(fi, . . . , f r ) and 



(33) 



/0 1 

1 



\ x ai a2 

Proof: Let v = (b, a, 0, . . . , 0). Then 



\ 



1 

1 

a r _ 2 / 



(34) 



r-2 



(35) 



( c Cfi)> • • • ! c (/r)) = vQ(x + A) . 
Now fi'(x) = J B(x)fi(x). Moreover, for j < r-2, = (0, . . . , 0, b, a, 0, . . . 

Thus 



But 



( V 
vB 



Wr(c(f ,),..., c(f r )) 



( v \ 

vB(x + A) 



n(x + A) . 



(36) 



\vB(x + xy- 1 J 

bl + aB, so the result follows. 



□ 



\vB r ~ l J 
Now we give the 

Proof of Theorem \2.4\ : For i = 1, . . . , n and j = 0, ... n — 1, let Aij be 

the r x r matrix 

:= c^-J + A[S(x + AO , (37) 

where 

&ij ■= liK • ( 38 ) 

Given any / e ker(L ), 

(c(f),PMf)),---,PN(c(f))) = (C(/),...,C(/^))), 

c(z/),...,c(^)),..., 
c(z n f)) . (39) 
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±1 



It then follows from Lemma |2.5| that 
Ku 

where 



A U 
-A Pn{u) 



(40) 



U 



( u \ 

P 1 (u) 



(41) 



\Pn-i{u)J 

A := (- 7o A£ - n\ n Q ~\ —Xfi,0, ... ,0, ... , - ln -iK-i ~ ™K-l ■ ■ ■ , 0) 



and 



A :-- 



r-2 

/ A 10 A 20 
A u 

V^l.n-l ... 



r-2 



A n Q \ 



An n— i ) 



(42) 



(43) 



It follows immediately that K has polynomial coefficients. Moreover, if we 
let A' be the n x n submatrix of A in which x occurs, then 



= (n(A. 4 - a,) j (n(*+*o) 



(44) 



Letting A" be the complementary (r — l)n x (r — l)n submatrix, we have 

(45) 



i4 w i = in(*-*j 



Thus 



14 = I II(* - A i) I + GOO 



(46) 



which proves assertion 2. To prove 3, let K^j be the coefficient of Pi r +j(u) 

(47) 



in the determinant Equation |4(| Then 

i^ = $>W M (Lo) . 



Since deg K n $ = deg |v4| = n by Equation |46], we must show that 

deg (K ifi ) < n and (48) 
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deg (Kij) < n for j > . 
Labeling the rows of A as A , . . . A N _i, 



±\A 



(m) | 



(49) 



(50) 



where m = ir + j and At m ) is the matrix obtained from A by replacing 
A m with A. Then the same argument that proves assertion 2 proves that 
deg (i£j,o) < n. It also proves that deg {KiJ) < n if < j < r — 1. Indeed, 
in that case the submatrix of A which gave us the coefficient of x n has now 
acquired a repeated row. But if j = r — 1, then we have removed one of the 
rows in which x occurred, so in that case the degree also drops. □ 



2. 6 Examples 



It follows from the proof of Theorem 2.4 that one may easily compute K 



and . In the last column of determinant Equation f|0], replace P r j + k(u) 
with z^ h . The determinant is then a polynomial in x, z, £. To get K, 
we make the replacement z^ k i— > d k L-Q, and to get X b we make the same 
replacement after first switching x and z. 

Case n = 1 and r = 2: In this case, we have c = 5a ° (d + 7) and 
Ln = d 2 — x. The determinant is 



-£ + 7 + 7 2 A-A 2 -Ax-(7 2 -A)^ + xz (51) 



7 1 1 
A + x 7 £ 
1 + 7A A z 

This gives 

= ( x + A- 7 2 )a 2 -9+(7 2 -A-x)x + 7 + 7 2 A-A 2 -Ax ,(52) 
K b = ( x - A)5 2 - d + (A - x)x + 7 + 7 2 A - A 2 - ( 7 2 - X)x . (53) 

The point, which we will generalize, is that is obtained from K by the 
involutive map 

(A, 7 )^(7 2 -A, 7 ) • 
In other words, K c = K c p for some 

Case n — 1 and r = 3: When r = 3, the vacuum has a parameter, 

L n = <9 3 — ad — a; . 



(54) 



-0 

With c as in the previous example, the relevant determinant is 
7 10 1 

a-7 2 + a7A-7 3 A-A 2 -Ax + 



(55) 



7 
A + x a 
1 + 7A A 



1 e 

7 e 2 

z 



7C-^ 2 -«7-z + 7 3 z + Az + a;2; (56) 
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This gives 

K = (x - 07 + 7 3 + A)<9 3 - d 2 + (7 + a 2 7 - cry 3 - aX - ax)d 
+a — 7 2 + cry A — 7 3 A — A 2 + cryx — j 3 x — 2Xx — x 2 
= (x-A)9 3 -9 2 + (7 + aA-ax)9 + a- 7 2 + a 7 A- 

7 3 A-A 2 -a7x + 7 3 a: + 2Az-x 2 (57) 

These correspond under the involution 

(A, 7 )^(a 7 -7 3 -A,7) • (58) 

3 Bispectrality 



Our aim now is to prove that for K = Kq as defined in Equation the 
rings Rr,l an d Rr\l are nontrivial. First consider the stablizer of C in 
the polynomial ring 

R c ■= {p( z ) g C[z}\c o p e C for all cG C} . (59) 

It follows from the formula 

5\°d z o p(z) = p(X) oS x od z + p(X) o 5\ (60) 

that 

R c = {p(z)\p'(X l )=0, l<i<n}. (61) 
Lemma 3.1: Rr,l — Rc- 

Proof: It is a standard fact that the left ideal T>K is precisely the set 
of operators which annihilate the kernel of K. In particular, p(z) belongs 
to Rr,l if an d only if Kp(L )(c(f)) = for every / G ker(L ) and c G 
C. Since Kp(L )(c(f)) = Kc(p(z)f), this is the same has saying C o 



p(z)(ker(Lo)) C C(ker(Lo)). By Lemma 2T , this occurs if and only if 



C o p(z) cC. □ 



It now follows from (|6lf) that Rr,l C and in fact contains elements 
of every sufficiently high degree. Consequently, 

K RM := [Kp{L Q )K- l \p G R RM ) (62) 

is a commutative ring of ordinary differential operators of rank r. We 



remark that this ring is the iterated Darboux transformation JIZ/ of C [Lq] 
by the zero eigenfunctions cf> s . 
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In fact, it is well-known that for any point C G Gr^S), Rc is non- 
trivial. Thus, to prove that Rx b c l * s a ^ so non-trivial it suffices to find a 
space C 13 G Grh(S) such that Kq = K C fs. 
Definition 3.2: Let C 13 C S be the subspace 

:= {c G S\c{f) G keri? b V/ G ker L }. 

An important alternative definition of C 13 is given as follows. 

Lemma 3.3: C p = {c G S\c o K(z, d z ) = 0}. 

Proof: Using the definition (H) of the involution b, c belongs to C^ 3 if and 
only if 

= K'oc(f) 
= coK\f) 

= coK(z,d z ){f) . (63) 
By Lemma |2.1| , this occurs if and only if c o K(z, d z ) = 0. □ 
Lemma 3.4: C 13 is at most n- dimensional. 



Proof: By Theorem |2]4] (3), K b has degree jV = rn. Thus, its kernel 
is at most iV dimensional. However, by Proposition |2.3| , the vector space 
{c(f)\c G C 13 f G ker L } has dimension r ■ dimC 3 . Since this is contained 
in the kernel of by definition, dimC 13 < n. □ 

Now we make a general observation about Wronskians, which will en- 
able us to conclude that the dimension of C 13 is exactly n. 

Theorem 3.5: Let g := (gi(z), . . . ,g m (z)) be a vector of analytic functions 
of z. Let w(z) := \Wr(gi(z), . . . ,g m (z))\ and let Q(z, d z ) be the differential 
operator whose application to an arbitrary function u(z) is given by Qu = 
\Wr(gi(z), . . . ,g m (z),u)\. Given A G C such that w(X) = and w'(X) ^ ; 
there exists a unique distribution of the form c = 5\ o (<9 + 7) such that 
c o Q = 0. Moreover, if we set 



Q = w {z)d™ - w\z)d^- 1 + v{z)8™- 2 + 0(a™- 3 ) , (64) 

then 

w'(X) 



v(X)-w"(X) 

7 = 77T^ • ( 65 ) 



Proof: We have 

H u 

(m+l), 



A x :=(doQ) (u) 



z=\ 



g(m + l) u( m+1 >(X) 



(66) 
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and 

A 2 := Q{u) 
where g(j) is the 1 x m vector 



H u 

g{m) u^ m \\) 



g(J) '■= T7^' • • • >9*> 



(67) 



(61 



H is the m x m matrix 



H :-- 




(69) 



and u is the m x 1 vector 



it :- 



/ «(A) \ 
tt'(A) 



(m-l) 



(70) 



(A)/ 



Since \H\ = 0, the "bottom right" corner of each matrix above does not 
affect the determinant and 



Ax + 7A2 



H u 
g(m + 1) + ^g{m) 



(71) 



This determinant is zero (for arbitrary u) if and only if g(m+l)+jg(m) G V 
for 

V := (g(s)\0 < s < m- 1) C C m . 



Furthermore, since w'(X) ^ 0, the vectors g(s) for s = 0, 
C m . Consequently, there is a unique 7 such that 



m-2 



<f( OT + 1) = -jg{m) + a i9(i) 

i=0 



(72) 
m span all of 

(73) 



Such a 7 gives us our c. On the other hand, a direct calculation shows that 
the highest order terms of c\ o Q vanish only if 7 has the form (|65D. □ 

In the rank one case, the bispectral involution was shown to exchange 
the roles of the r-function, which is zero when the x-orbit leaves the big cell 
in the grassmannian, and the polynomial q(z) which is zero at the singular 
points of the spectral curve \T3rf . This pair of polynomials will be seen 
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to satisfy the same relationship and will also play an important role in the 
present paper. 

Definition 3.6: Let r{x) = tc(x) be the coefficient of the d N in K. In 
particular, 

t(x) := -, — ^\A\ (74) 

is a monic polynomial of degree n where A is as defined by ( |37"D and fl43|) . 
Let q c (z) = q(z) :— Ui=i(z — A*). 

Let D(n) be the subset of D consisting of subspaces of dimension n. 
We immediately have 

Corollary 3.7: If C G D(n) and r has distinct roots, then C 13 G D(n). 
Moreover, C 13 has support at the roots of the polynomial r. 

Now we can deduce the main results. For convenience, we denote K 13 : = 
K c i3, t^(x) = Tq0 and denote by D the subset 

D := {C G D|r c has distinct roots} . (75) 

Theorem 3.8: Assume C G D . Then 

1. t 13 = q 

2. i? b = K 13 . 

Consequently, one may conclude that f3 is an involution on D . 

Proof: By definition, belongs to the ideal of operators annihilating 
c(f) for / G ker(Lo) and c G C 13 . By Corollary |3.7| and Theorem |2.4] , both 
and K 13 have degree N. Thus there exists a non-zero rational function 
/i £ C(i) such that 

K b = fxRf 3 . (76) 

Observe next that that the the coefficients of K b have degree at most n. 
This follows from the definition of b and the fact that K has order N. 
Moreover, c o K b = for all c G C. It follows that the leading coefficient 
of Sy vanishes on the support of C. This coefficient is easily seen to be 
monic, whence 

K b = q(x)8 N + b(x)d N - 1 + 0(d N ~ 2 ) . (77) 
for some polynomial b(x) of degree at most n. We also have 

K 13 = r p {x)d N - r' p {x)d N ^ + 0(d N - 2 ) . (78) 
Thus both claims are proved once we show that fj, — 1. 
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Again since c o K b = 0, 

6(A) = -q'(X) (79) 

for all A in the support of C. Thus, there exists a constant ( G C, such 
that 

b(x) = -q'(x) + £q(x) . (80) 



Then 



pr p = q 



f*r* = q'-Hq- (81) 

Taking quotients and integrating, we find that there exists p G C* such 
that 

t P {x) = pq(x)e^ x . (82) 

Since r 13 and q are polynomials, p — 1 and £ = 0. 

For any c G C, c(f) G ker K by definition of X. Furthermore, since 
= g, it has distinct roots implying that C G Do- Therefore, by the 
results above, (K^f = (K b f = K and so c(f) G ker (K /3 f. Thus C C 
{C^Y . Since both spaces have dimension n, they are equal. □ 

Theorem 3.9: For any C G D and K (= Kq ), the ring IZr l ^ s a bispectral 
ring of operators. 

Proof: By Proposition [L2| it is sufficient to show that Rr,l an d Rr^> l 



both contain non-constant polynomials. It follows from Lemma [3.1| that 
Rk,l is non-trivial. One can show that if Qt is a family of elements in the 
Weyl algebra depending continuously on a parameter t, Q t having a fixed 
order in both d and x for all t, then the property of RQ t) L being nontrivial 
is preserved under limits with respect to t. Thus we reduce to the case that 



r has distinct roots. In such a case, Theorem 3.8 shows that K = Kqp 



for C 13 G D and so, once again, Lemma |3J] demonstrates that Rr\l i s 



non-trivial. □ 
Definition 3.10: For any / G kerL , let 

i c (x,z):=-^-Kf(x,z). (83) 

The map / i— > f c takes ker L to an r-dimensional space of common 
eigenfunctions for the operators Q p = Kp{L§)K~ x G TZr Lo , satisfying 



Q p f c (x,z) = p(z)f c (x,z). 



(84) 
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As with the rank one bispectral involution [JIJ, the action of (3 exchanges 
the roles of the spectral and spatial parameter in the eigenfunctions, as 
demonstrated by the following result. 

Proposition 3.11: For any f G kerL and any C G D the eigenfunctions 
fc and ic0 are related by the formula 

fc{x,z) = i c p{z,x). (85) 



Proof: 



fc(x,z) = -}—Kf(x,z) 
q(z) 



' -Kf(x,z) 



r(x)q(z) 

1 K b (z,d z )f(x,z) 
r(x)q(z) 

1 q(z)K^z,d 2 )f(x,z) 
T{x)q{z) 

1 



qf>(x) 



K^z,d z )f(x,z) 



fcp{z,x). (86) 

□ 



4 True Rank 

The ring of operators TZ^ Lo clearly has rank r. However, if this ring is 
contained in a larger commutative ring of lower rank, it is said to only 
have "fake" rank r [7^ /. For example, if some other operator of order 
relatively prime to r commutes with the elements of TZ^ L , then it actually 
has "true" rank one. Since the bispectral operators contained in rank one 
bispectral rings have already been identified by Wilson /[73|/, it is useful to 
note that the ring 7Zk,l ^ s n °t contained in a commutative ring of rank 
smaller than r. 

Definition 4.1: The true rank of an ordinary differential operator is the 
rank of its centralizer in the ring of all ordinary differential operators. The 
true rank of a commutative ring of operators is the true rank of any of its 
elements. 



Lemma 4.2: The centralizer of L in the ring of ordinary differential oper- 
ators is the polynomial ring C[L ]. In particular, L has true rank r. 
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Proof: Since L is contained in the Weyl algebra of ordinary differential 
operators having polynomial coefficients, its centralizer is as well. Since b 
is an algebra antiautomorphism, it suffices to prove the result for Lq . But 
Lq = x, and it is well-known that the centralizer of x in the Weyl algebra 
is C[x]. □ 

Lemma 4.3: If X = Y X Y 2 and X = Y 2 Y\ then X and X have the same true 
rank. 

Proof: If Q is an operator commuting with X, then 

= Y 2 [Q,X]Y X 

= Y 2 QY 1 Y 2 Y 1 -Y 2 Y 1 Y 2 QY 1 

= [Y 2 QY X ,X]. (87) 



Let r be the true rank of X. Then, by (|87|), we have ord(Y 2 QYi) = mod r. 
But, ord{Y 2 QYi) = ord{Y 2 Yi) + ord(Q) and since ord(Y 2 Yi) = mod r we 
conclude that ord(Q) = mod r. Therefore, the true rank of X divides 
the true rank of X. Then, by symmetry, the true ranks are equal. □ 

Note: Lemma [4.3| generalizes a previous result demonstrating that 
Darboux transformations preserve true rank in the case r = 2. 

Proposition 4.4: The ring TZk^o °f ordinary differential operators has true 
rank r. 

Proof: The operator q 2 (Lo) G C[Lo] has true rank r, according to 
Lemma [4.2| . Furthermore, 

q 2 (L )c(f) = c(q 2 (z)f) =0 VceC, fe ker L . (88) 

Thus, in particular, q 2 (Lo) = QK for some Q G T>. Then the operator 
Q q 2 = Kq 2 (L Q )K^ 1 G TtK,L is given by KQ and has the same true rank 
as q 2 {L ). □ 



5 Discussion 

5. 1 Examples 

The following example demonstrates the construction of bispectral algebras 
of rank r = 2 in the case n = 1. For r = 2, there is only one choice for a 
generalized Airy vacuum: L = d 2 — x. Thus, we may choose the classical 
Airy functions f\(x) = Ai(x) and f 2 (x) = Bi(x) as a basis for kerLo- 
Consider the one-dimensional space C = (5q o (<9 2 + 7)). A polynomial p(z) 
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stabilizes c if and only if p'(0) = 0. Consequently, Rc = C[z 2 ,z 3 ]. From 
Equation fl52|) r(x) = x — 7 2 and 



7^ — x 



7^ — x 



Then from Definition 3.10 



fc{x,z) = (1 
and from Equation 



7 



z(7 2 — x) 



)f{x,z) + 



z{^ 2 — x) 



59) 



/'(x 5 z) (90) 



(91) 



where 



9 4 + 



2 



2 (-2 - 7 4 x + 27V - x 3 ) 
(7 2 - x) 2 

-2 (-4 + 2 7 3 - 7 6 - 2 7 x + 37 4 x - 37V + x 3 ^ 



+x^ 



(-7 2 + x) 3 
47 



-a 



+ 



(-7 2 + x) 4 (- 7 2 + x) d -7 J + x 



and 

u 



d b + 



3 (-2 - 7 4 x + 27 2 x 2 - x 3N| 



+ ^-6 + 
+ ( 3x 



(7 2 — xY 
24 



-d 4 



+ 



6",' 



(-7 2 + x) 3 (-7 2 + ^) 2 , 
72 I87 



d 3 



+ 



6 7 2 



+ 



(- 7 2 +x) 4 (-7 2 + x) d (-7 2 +x) z -7 2 + ^ 



,9 2 



/ 144 
+ 6z + nr + 



367 



12 7 2 



3 (3 + 2 7 3 ) 67 



(-7 2 + x) 5 (- 7 2 + x) 4 (- 7 2 + x) 3 (- 7 2 + x) 2 -7 2 + 



-1 - x J 



144 



367 



+ 



12 7 2 



(— 7 2 + x) (— 7 2 + x) (— 7 2 + x) 



3 (3 + 2 7 3 ) 



3-) 



37 2 



(-7 2 + x) 3 (- 7 2 + x) 2 -7 2 + * 

ion (Q) we se 
i? c , = C[(z-7 2 ) 2 ^-7 2 ) 3 ]and 



From Equation (|53| ) we see that = (6^2 o (c^ + 7)) and so t^(x) = x, 



1 

— i 

x 



K p = a 2 _ l 9 - x + 2 _ 7 2 



X 



(92) 
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Most significantly, computing ^(cc, z) := ^^/(x, z) it is easily deter- 
mined that 

fci>(x,z) = ?c{z,x) (93) 

and so the operators in 1Z c /3 := KC[(L — 7 2 ) 2 , (L — 7 2 ) 3 ]if _1 demonstrate 
the bispectrality of TZ^ Lq . 

In the case that 7 = 0, C and C 13 are the same and so this is a fixed 
point of the involution. This example is discussed by Griinbaum as 
a solution of the KP hierarchy. For arbitrary 7, the time dependent KP 
solution corresponding to this ring is a non- vanishing rational solution jj^J 
and is described in j^j . 

5.2 Sato Grassmannian and KP Flows 

The composition of the Wilson's involution on Gr 1 with the KP flows was 
shown to be a non-isospectral symmetry of the KP hierarchy and a lin- 
earizing map of the Calogero-Moser Particle System Qj. The higher rank 
involution presented here may similarly have interesting dynamic proper- 
ties. Therefore, the embedding of Grh{S) into the Sato grassmannian Gr r 
used in constructions of rank r KP solutions ]7I}/ may be of interest. 

Essentially, the construction above is a special case of the "dual con- 
struction" j^j which associates a point in Gr r to subspaces of S r whose 
dimension is divisible by r. Let 

F(x) :=(Wr(f 1 (x),...,f r (x)))- 1 

and ^Fj{x) denote the j column of this matrix. Define the composition of 
an r-vector v = (vi(z), . . . , v r (z)) with an element c G S to be the element 
co v = (c o vi, . . . , c o v r ) G S r : 

(c o v) (w) := c(v • w). (94) 

Then C C S T is the iV dimensional subspace 

C := (co^ j (z)\ce C 1 < j <r). (95) 

To the vector space C, we associate the pointlll W G Gr r 

W := -}-rVc (96) 



1 In the case that q(z) has a root on the unit circle S 1 , the point W achieved in this 
way is not contained in Gr r as formulated in JT^/ since its elements were taken to be 
L 2 (S' 1 ). This subtlety need not concern us here, however. 
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where Vg is the null space of C in (C[z]) r and the overline indicates Hilbert 
closure in H r = L 2 (S l ) r . 

Proposition 5.3: The vector Baker function ip^(x,z) and the eigenfunc- 
tions functions fj(x,z) := -^K fj(x, z) corresponding to C are related by 
the formula 

rfarfaz) = (fi(x, z),. . . , f r (x, z)) ■ F^ 1 (z). (97) 

Furthermore, the ring of ordinary differential operators associated to W is 
the ring Urm . 

Then, if we denote by W the point in Gr r corresponding to C and by 
W 13 the point corresponding to C 13 , Proposition |3.11| implies 

Corollary 5.4: The action of the involution (3 on vector baker functions is 
given by the formula 

ihy(x, z) = ijj w0 (z, x) ■ F-\x) ■ F(z). (98) 



5. 5 Bispectral Flows 

In light of the remarks of the preceding section, it follows from the results of 
/H/ that composing the distributions in C with the function e tz induces 



z k 

the rk th flow of the KP hierarchy. In particular, C = K(t)L\' 7 K~ l {t) 
satisfies the equation 

d 



For C E D and t sufficiently small, C o e tz G D . Consequently, one 
may conclude that the KP flows whose indices are zero modulo r are a local 
symmetry of these bispectral solutions. Furthermore, as in jffl, the com- 
position of the KP flows with the bispectral involution is a non-isospectral 
symmetry. Finally, it can be shown that the bispectral involution is a sym- 
plectic map linearizing the Hamiltonian system determining the motion of 
the poles of these KP solutions. (We intend to discuss these results further 
in an upcoming paper.) 

5.6 General Remarks 

Associated to any commutative ring of ordinary differential operators of 
rank r is its geometric spectral data, which are a rank r vector bundle over 
a complex projective curve J| ; [71|/ . The geometric spectral data corre- 
sponding to rank one bispectral rings are "unicursal" rational curves with 
any line bundle fiTW ■ That is, any rational curve having singularities in 
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the form of cusps with any line bundle is the geometric spectral data of 
some bispectral ring of operators. The spectral curves associated to the 
bispectral rings constructed in this paper are still rational curves with sin- 
gularities only in the form of cusps. However, the bispectral rings presented 
here do not correspond to arbitrary rank r vector bundles, but rather only 
very special bundles. Indeed, for given n and r, our construction depends 
on n + r — 2 parameters, which is in general fewer than the number of 
moduli for a bundle of rank r on a singular rational curve of genus n. 

Previous work has indicated a relationship between bispectrality and 
polynomial r-functions [7||/. The results above continue to support 
such a relationship (Definition |3.6|) . Furthermore, the bispectral involution 
involves not only an exchange of the variables x and z, but of the roles of 
the polynomials t(x) and q(z). Since the roots of r(x) indicate places that 
the orbit of the point W G Gr r leave the big cell, and the roots of q(z) 
indicate the singular points of the spectral curve, the bispectral involution 
can be roughly seen as an exchange of spectral curve with the x-orbit. 

The rings 1Z constructed in this paper contain operators normalized so 
that they are monic and have zero second coefficient. A change in variable 
or conjugation by a non-constant function would change this normalization 
but preserve bispectrality jl^/. 

It has been noted previously ^ that in the case of the known bispec- 
tral operators associated to the KP hierarchy, one is able to determine a 
corresponding operator in z whose eigenvalue is any antiderivative of the 
r function. The same can be shown to be true in the present case. In par- 
ticular, for every p G Rcp, there is an operator Q p (z, d z ) for which fc( x , z ) 
is an eigenfunction with eigenvalue p(x). However, since q l3 {z) = t(z), we 



have by Equation |61] that R c p is exactly the ring of polynomials whose 
derivatives are divisible by r. 

Finally, in light of Wilson's result, it is natural to conjecture that given 
L , one may associate a bispectral ring to any finite dimensional homo- 
geneous subspace C C S by defining Hk,l as above. In particular, we 



conjecture that if we define as in Definition EO for any C G Gr^^S) 



there always exists a nonzero \x G C(x) such that 

i? b = fi(x)K p , (99) 

which would imply that 

r k\l = Rkp,l = Re? C . (100) 



For instance, suppose r = 2,C = (5\ 1 o (<9 + 71), £\ 2 (9 + 72)), Ai 7^ A 2 , 
but C does not belong to D . Then one may check directly that there 
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are two possibilities for C 13 . Letting A be the unique (repeated) root of r, 
C 13 is spanned by two distributions supported at A, of orders either (3, 0) 
or (2,1). In both cases, K@ has order 4, as does K < ". This implies that 
Equation d99|) holds. 

Note: After the completion of this work, a preprint /fj/ was made available 
which presents similar results from a Lie theoretic point of view. 
Acknowledgements The authors are grateful for friendly and helpful 
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